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Abstract 

We further analyse, at next-to-leading logs level, the form of fc-factorization and 
the definition of impact factors previously proposed by one of us, and we generalize 
them to the case of hard colourless probes. We then calculate the finite one-loop 
corrections to quark and gluon impact factors and we find them universal, and given 
by the same /C factor which occurs in the soft timelike splitting functions. 
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1 Introduction 



The rise of small- a; structure functions at HERA ^ has emphasized the need of under- 
standing high-energy QCD 0] at next-to- leading logx (NLx) level |13|-||1CI||. 

The general purpose of such calculations is to derive the NLx corrections to both the 



BFKL equation and to the resummed anomalous dimensions JTOl in hard processes 



with one or more hard scales. Therefore, we need to compute on one hand the NLx 



BFKL kernel [O, O, which is supposed to be independent of the external probes, and 



on the other hand the impact factors [|rT[|, which characterize the probe. 

The separation of impact factors and kernel is made on the basis of high-energy k- 
dependent factorization (Sec. (|])), or in short fc-factorization [Q, which is therefore to be 
extended at NLx level. 

The purpose of the present paper is to further analyse a form of fc-factorization pro- 



posed by one of us ||TT[, by defining the impact factors in a proper factorization scheme 



and by computing the one-loop corrections to them in the case of partonic probes. 

We have emphasized the use of a proper factorization scheme, because the separation 
of the cross section in impact factor contributions and gluon Green's function suffers 
of some ambiguity, analogous to the one between coefficient functions and anomalous 
dimensions in collinear factorization. In fact, the subtraction of the leading logs terms 
involves a prescription which is not unique, not only for choosing the scale of s, but also 
for the form of Q^^ at finite energies. Therefore, some probe-independent NLx terms can 
be attributed to either the impact factors or to the Green's function, depending on the 
scheme being adopted. 

A definition of impact factors on the basis of a proper scale choice for subtracting the 
leading term was given [|ll|] for the case of dijet production, and was later used ||13[ in 
order to extract the NLx kernel []T2|, |TB[. Here we further discuss this definition and the 
ensuing factorization scheme in Sec. (|^). First, we emphasize the requirements that the 
subtracted leading term should satisfy at finite energies for the impact factors to be well 
defined (e.g., without spurious infrared divergences). 
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Furthermore, we generalize the definition above to the case of colourless probes, which 
is the most important one for applications to DIS and to heavy quark processes. We argue 
on this basis that the NLx kernel extracted before in a partonic process applies to 



physical hard processes as well, in the given factorization scheme. 

We then proceed to the actual calculation of partonic impact factors in Sees. (|,^. 
The latter still contain some collinear divergences due to the massless initial partons, 
which are those not subtracted out with the leading terms. We then factor them out in 
proper splitting functions, and we end up with some finite one-loop corrections not only 
for the quark |Tl|, but also for the gluon initial state. 



It is amusing to note that such finite renormalization factors (Sec. (^) are universal, 
i.e. are the same for quark, antiquark and gluon, and are given in terms of the coefficient 
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in the fc-factorization scheme of Ref. [ITTI. In other schemes 



such finite renormalizations are different, but still independent of the probe. 

A similar property, and the same coefficient /C, were noticed to occur in the one-loop 
corrections to the 1/(1 — z) part of the timelike DGLAP splitting functions. In this 
case the universality feature is probably to be ascribed to soft gluon properties [T^ which 



appear to determine not only single parton emission, but correlations as well. Apparently, 
a related argument can be applied to gluon-Regge exchanges, along the lines suggested 



by Ref. [15 



We finally discuss our results in Sec. (0), by giving some calculational details in Ap- 
pendices A and B. 



2 k-Factorization in dijet production 

Most of the NLx calculations needed for the definition of the NLx kernel have been 
performed 0, 1^ in the case of high-energy parton-parton scattering, in which no physical 
hard scale is present. While the total cross section of this process has severe power-like 
Coulomb singularities, it is hoped that by fixing the transverse momenta ki , k2 of the 



fragmentation jets (corresponding to the virtualities of the exchanged gluons, Fig. (|l]), 
one is able to define a two scale hard process. A similar but not identical procedure was 
devised by Mueller and Navelet [1^ in hadron-hadron scattering (cfr. Sec. (P)). 

We then consider high-energy scattering of two partons a,b (a,b=q (quark), g(gluon)) 
with momenta pi,P2 and we factorize it in impact factors and gluon Green's function 
contributions as in Fig. (||), with the definition of Fig. (H). 




Figure 1: Double k-factorization of dijet differential cross section. 



By definition, the impact factors are free of high-energy gluon exchanges, which are 
subtracted out as in Fig. (0), but still have collinear singularities due to the initial mass- 
less partons which need to be factored out. Therefore, the Regge-gluon exchanges are 
incorporated in the Green's function Q^^ of Fig. (|l]), and define in a gauge- invariant way 
the virtualities ki and ^2 which play the role of hard scales of the process. 

Formally, we can write [|ll| the colour averaged differential cross section 

where we have introduced, besides the /I's and Q^^, also the energy-scale so(fci, ^2) which 
is normalized by the condition so{ki,ki) = k1 and will be mostly chosen to be sq = 



^1^2 (ki = \ki\), for the reasons explained below |I2| , 0. We adopt d[A;] = d k/n^^ as 
transverse space measure. The validity of Eq. (|2.1| ) at NLx level has been so far checked 
at one-loop and two-loop |13| orders and needs be established at higher orders. 
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Figure 2: Diagrammatic representation of the impact factor. In the sum over the possible final 
states, the outgoing parton index c denotes the gluon g for a = q while it denotes the quark q 
or the softer gluon for a = g. 



Let us start reviewing the leading logx (Lx) results and the NLx ones of Refs. fTT 



and 13 



We use the notation pi,P2 (Ai, A2) for the initial parton's momenta (helicities) and the 
indices 3,4 (possibly 5) for the final ones, with the Sudakov parametrization 

ki=Pi-P3 = zipi - ^—-P2 + kii. , 

1 - zi)s 



k2= P2- Pi 



-Pl + Z2P2 + /C2± 



(1 - Z2)S' 

where we have introduced Sudakov variables Zi and transverse spacelike vectors fcjj, per- 
pendicular to the plane {pi,P2) '■ Pi-ki± = = p2-ki± with D — 2 euclidean components 



At Born level, the cross section ( p^.l] ) is just non-abelian Coulomb scattering (Fig. ^, 
corresponding to the high-energy amplitude 



(2.2) 



where helicity conserving factors are exhibited, but colour indices are understood. 

Upon helicity and colour averaging and by using the (4 + 2e)-dimensional phase space 
for two particle final states 

TT 1 



2) 



2+e 



(2.3) 



Figure 3: Born diagram of high energy scattering. 



Eq. ( |2.2|) yields the differential cross section 



d[k] - 1 



C^aa(l^) =h^^\k)h^^\k) , (2.4) 



-1 4 UnY^^ 



2N^ 3 ' ^ ' ' r(l-£ 

in terms of the dimensionless coupling constant 



rv 

(47r)i+- 

and of the Born impact factors 

being the renormalization scale. 

At one-loop order, real emission and virtual contributions must be added to the Born 
cross section. At Lx level, the high-energy kinematics favours the emission of an additional 
gluon (Fig. ^) of momentum P5 = q and polarization e in the central region ^ q~ ^ \q\ ■ 
The corresponding amplitude is 

■^il=9' ^f^J'^^^, t:tir'''J-^KAx. (2.6) 
where q^ = + ki^ and 

JM(A;,, k2) = -A;f + k^ + ^(g2 _ kj) - ^(g2 _ kl) (2.7) 

Pi-q P2-q 



is the Lipatov current [Q. 



Figure 4: Real gluon emission in central region. 



After polarization and colour averaging, and by using the three-body phase space 



^(3) 



TT dzi d[ki] d[k2\ 



(2.8) 



we obtain the differential cross section 



real 



which, upon the Zi integration with infrared boundary Zi > Sq/s, becomes 

-/.(°)(fco/^r(fe2) , log- . (2.10) 



d[fci]d[fc2] / a V X, b V ^^g2r(l_^)^2. 

/ real 

At this stage, the scale Sq is totally undetermined (pure phase space would yield 
^ < zi < 1 — J-) because the NLx constant has not been calculated yet (cfr. SecH). 

The Lx virtual one-loop corrections give rise to reggeization of the exchanged gluons, 
so that their propagator is modified by a factor 

0^"^ = log ^ + Oia',) (2.11) 

where 

^ ^1^)- (4vr)2+e'' 7 p2(fc_p)2 - 2er(l + 2e)Uv ' 
Therefore, the one-loop contribution to the elastic cross-section becomes 



d[k] / d[fc] ' ' ^fc 

virtual 



2cu«(fc2)log^ (2.13) 



and the total leading contribution, up to one-loop order, becomes 



da 



ab 



d[fci] d[k2] 



S[q] + asKo{ki,k2)\og- 



So 



(2.14) 



where 



asKo{ki,k2 



a. 



g2r(l -£)/i2^ 

is the leading BFKL kernel B. 



+ 2u^'\kl)6[q] ■ 6[q] = n'+^6'+'%q) (2.15) 



The result ( |2.14| ) checks with the representation ( |2.1| ) by defining the Green's function 



(2.16) 



which satisfies the Lx BFKL equation. Next, we shall consider NLx corrections. 



3 Impact factors at NLx accuracy 

At one-loop NLx level, we should calculate the constant contributions to the differential 
cross section (p.l|) , which require an accurate treatment of the fragmentation region. 
Constant terms arise from both the one-loop corrections to the impact factors and possible 
constant kernels in the gluon Green's function. In fact, the general NLx form of Q^^ is |T^ 

-1 



UJ 



[I + cIsHr) (3.1) 



where we have introduced the operator factors Hi, Hr and the NLx kernel 

K^^ = -haAog^K, + asK, ; = (3.2) 



which consists in the running coupling contribution, proportional to Kq with a /3-function 
coefficient, and in the renormalization scale independent part agKi. 

The one-loop expansion of Eq. (|2.1|) , taking into account the definition ( p.l|) , reads 
/i(°)(fci) asK,{k,M) h^"\k2) log 

so{ki, k2) 

+ as[h<^^\k,)5[q] + h<^^\k,)HLik,,k2))h^^\k2) + a^b . (3.3) 



We see, therefore, that the identification of the one- loop impact factors h^a \ h^^^ is affected 
by a factorization scheme ambiguity, due to both the choice of the scale Sq and of the 
kernel {Hr). 



In Ref. [O we have defined the impact factors by the following factorization procedure: 



{i) Subtract the leading term with a reference scale, consistent with the infrared and 
collinear properties of the process, for which we can set = Hji = 0; 

(u) Interpret the remaining constant at single-fc factorization level as the one-loop cor- 
rection h^^\ 

In order to perform the subtraction (z) it is convenient to define the fragmentation ver- 
texQFa(zi, ki, ^2) which contains all (real and virtual) one-loop squared matrix elements, 
after factorization of the exchanged Regge-gluon k2. In the real emission contribution 
to -Fa, we fix the momentum fraction Zi and the transverse momentum q of the emitted 
parton c. In order to avoid double counting we also set c = g for a = q and c = q or c = 



the softer gluon in the case a = g. Then in Ref. the leading term was subtracted out 



as follows (here and in the sequel q = \q\, = \ki\): 

FM, ku k,) - a./i(°)(fci) ^^^2r(7!t)^2^ = ^sFi'\z,, k,, k2) 
and the impact factor contribution was defined by zi and fci integration: 



(3.4) 



h^^\k2) = j d[k,]j\z,F^'\z,,k,,k2 



(3.5) 



Note that the theta-function in the subtraction term in Eq. ( p.4|) cuts off the leading 
contribution l/^iq^ for q < A^i, by the angular ordering constraint Zi < q/ki. This insures 
that Fi^'* be both integrable at 2:1 = and without spurious g = singularities. 

In fact, the fragmentation function is expected to have the leading 1/ Zi singularity 

for zi — > 0, but no = singularity^. This coherence effect, first noticed by one of us 

^In this paper Fa differs from that of Ref. |Q by a normahzation factor as/ ziq^ . 

^The l/<7^ terms, pertaining to the Sudakov fragmentation region, actually cancel out with virtual 



corrections in Eq. (3.5) (App. g) 



||17| , is to be inserted in the subtraction term also, in order to avoid spurious singularities 
in Eq. 

With the definition (3^), the impact factor (3^) is expected to have the = 



collinear singularities of only, except the ones subtracted out with the leading term. In 
the quark case, the left over collinear singularity is the one provided by the nonsingular 
part of the zero- moment Pgq{uj = 0) (see Sec. @). 

By integrating Eq. (|3.4|) in the fragmentation region of a (-^ < 2i < 1) we obtain 



dziF^{zi,ki,k2) = as . 

q/^s q^r(l - 



^ Max(g, /ei) ~'~ ^ 



(3.6) 



where /la, which is characterized by its vanishing at g = [|lT|, is related to the proper 
anomalous dimension for ki = 0, as shown in more explicit form in Sec. (^). Therefore, 
the definition ( ^.51 ) corresponds to setting = Hr = for the reference scale 

s> = Max(g, A;i)Max(g, /^a) (3.7) 

which takes into account the coherence effect mentioned before. 

A similar definition of the impact factors can be adopted for colourless sources also, as 
in double-DIS or quarkonium processes. In this case, there is an additional dependence 
on the u variable and on the external hard scale (s) so that, at leading level, the impact 
factors 

.2 ■ 



(3.i 



yield a nontrivial fci-dependence, which has been explicitly computed in several processes 
|1T8| . In Eq. (|3.8| ) the function /a has the role of setting ki of order Qi in the total cross 
section, and has no fc^ = singularity at all, because there is no initial colour charge^. 

Therefore, by translating angular ordering into the a;-space threshold factor (g/Zci)'^, 
the analogue of Eq. (|3l^) becomes 

F.Mi, fci, fcs) - ash'^^hQi, fci) (^] = a.Fi,l](Qi, fci, k^) (3.9) 



^We are indebted to V. Fadin for a discussion on this point. 
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with the definition 



h^^l{Qi,k2) = jd[k^]j\z,F^^{Q^,z,,k,,k2) . (3.10) 

In this case, no g = nor fci = singularities are expected in F3, so that it is again 
important that the leading term be subtracted out with the angular ordering constraint. 



We thus conclude that the definition of the impact factors in Eqs. ( |3.5|) and ( |3.10|) , on 



the basis of the subtraction in Eqs. ( |3.4| ) and ( ^.91 ) defines a self-consistent fc-factorization 
scheme of Eq. (|2.1|) for both coloured and colourless sources. 



The question then arises: since this factorization scheme singles out the scale s> of 
Eq. ( p.7|) , does this mean that Hi = Hr = in general? 



The answer to the question above depends on the identification of the scale sq, possibly 
different from s> , for which the representation ( p.l|) is assumed to be valid to all orders. 
Such all order identification has not been rigorously provided yet. However the scale s> 
has the defect of not being factorized in its ki,k2 dependence. This fact is argued to 
contradict multi-Regge factorization of production amplitudes |]19|, which implies short- 



range correlations of the fc's. Furthermore, the use of 0(2, 1) variables implies that 
Regge behaviour in the energy s should involve the boost 

co*(C) = 5^ . (3.11) 



thus suggesting the choice Sq = /ci/c2, adopted in Refs. [jT2| and ||13| 



If So = kik2 is assumed, then Eq. ( |3.6| ) implies that Hji = Uj^ = if is nonvanishing, 
and is given by the expression |]13| 



H{ki,k2) = ^logfe.fc, , iso = hk2) . (3.12) 

q H Ki 

By replacing Eq. ( p.l2|) in the form (|3.1|) of the gluon Green's function, the represen- 
tation ( |2.1| ) can be expanded at two-loop level in order to derive the NLx kernel Ki 



for the scale Sq = kik2, from known two-loop calculations |10]. The scale changes from 
So = kik2 to other factorized combinations, like Sq = kf {sq = A;|), which are relevant for 
kf > A;| {kl ^ kj), were described in Ref. ]T3[ . 
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An important point concerns the possible choice of different fc-factorization schemes. 
In fact the precise form of the subtraction term in Eqs. ( |3.4| ) and ( p.9| ) is not uniquely 
determined, away from zi = 0. Alternative forms exist, which reduce to the leading 
term for zi 0, and vanish for g — 0. Switching to such alternative forms defines an 
alternative A;-factorization scheme which differs from the above by a universal and finite 
impact factor renormalization. For instance, one can switch |Tl| from the reference scale 



s> to the scale Max(fc;^, fcg) by the change in the constant piece in Eq. 

6h{q, fci) = log ^ <dk2ki - log ^ Ogfci (3.13) 
which, upon fci-integration, implies an a-independent change in the impact factors 

5/.« = a.^^'(l)/i(°) . (3.14) 

Correspondingly, the gluon Green's function ( p.lD changes in the constant kernels by the 
quantities 

6HL{ki, fcs) = -6h{q, fci) , 6HR{ki, k^) = -6h{q, k^) (3.15) 
so as to keep the overall cross section unchanged. 



In the following, we apply the fc-factorization scheme defined in Eqs. (|3.4| ) and (|3.5| 
to the calculation of the one-loop corrections to the partonic impact factors. 



4 Quark-parton collision 



The calculation for the quark-parton scattering was done in Ref. [|T^ and here we sketch 
the main results. The real emission amplitude, assuming the gluon in the quark a frag- 
mentation region and adopting the Feynman gauge, is 

A^qgb = e^u,A%m ^g2pl^tl5x,x, (4.1) 
where ^4°^ is diagrammatically represented in Fig. (|^). 
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Figure 5: Real gluon emission amplitude in quark a fragmentation region. 



fi, a 



u, b 




+ 




+ 




Figure 6: Quark to gluon fragmentation tensor A"^. 



After a simple calculation, the differential cross section in this kinematical region 
results to be 



dcr, 



qgb 



dZid[fci]d[fc2] 



X 



f CAas (1 - zi)q-{q - Zifc a) , CpOs zjkl 



+ 



TT q^{q- Zik2Y 



(4.2) 



where 



'y,^{z,,e) = ^^[l + {l-z^f + ezl] 



(4.3) 



is related to the quark to gluon splitting function. We can easily check that dfXqgb matches 
dcTqg^ ( p.9|) in the central region zi <^ 1, so that we can take it as the right differential 
cross section in the whole positive rapidity range 



l/5 = log^^0 



Zl > 



(4.4) 
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In the remaining half phase-space < the cross section has the same expression pro- 
vided we exchange ki ^ k2 , zi ^ Z2, where Z2 is the momentum fraction of q with 
respect to p2- 



It is worth noting that the splitting function ( [4 .31 ) is factored out in Eq. ( |4.2| ) even 
outside the collinear regions ^ fc^ ~ k^ and fc^ ^ qr^ ^ k^ thus suggesting a smooth 



extrapolation between collinear and Regge regions |17, 21 



The correction to the cross section due to virtual emission, including subleading effects, 
for general parton-parton scattering, can be extracted from the amplitude of Ref. ^ 

s 



-Mab = 2.(^?t:5A3A.)(l + rW)l 

where, in the quark case. 



l+Lu{-t) log- 



-t 



1 + r(+) 



bb 



b"A4A2 



12 



36 



4 



N, V6 



Cf fl 3 ^ 
i^U-2+^^ 



(4.5) 



(4.6) 



is the quark-quark-reggeon vertex correction in the helicity conserving channel. Therefore, 
we have to add the virtual term h''q\ki) 2rq^''5(l — 5;i)5[qr] to the fragmentation vertex 

Fa(zi,fci,fc2). 

The Cp term of real emission is not singular at = and doesn't contribute to the 
logarithmic part of the cross section; therefore, it can be immediately integrated over zi 
and ki at fixed k2. Since the result exactly cancels the Cp term of virtual emission, we 
can write the total fragmentation vertex for the Ca term only: 

'11 



Fq(^i, fci, k2) = /if (fci) <! 6[q]6{l - z^)2u;^'\kl 

'?gq{zi,e) as (1 - 



+ 



12 



85 7r^\ n/ /I _ 5 

36^Ty^~iVeV6~18^ 
' 2:1^2 



r(l - e)iJ^^ q2 



(q - zik2 



(4.7) 



We can see how the features anticipated in Eqs. ( |3.4|) and ( p.5|) are realized by the 
explicit expression (|4.7|) . In fact, Fq has no q^ singularity, except at = 0, where it 
reduces to the leading term to be subtracted out. 



By using the procedure of Eqs. ( p.4|) -( pT6D we arrive at the definition 

' ^ hq{qM)= / dziF«(^i,fci,fc2) (4.^ 

13 



q^T{l - £)/i2^ 



in which hc^{q, ki) vanishes at q = 0. By integration at fixed k2, we obtain (App. |BD the 
one-loop correction to the quark impact factor, as follows: 



^2 



/c 



67 
18 



TT 

"6 



11 
6" ' 

9 



3Nr 



+ 



1 \ fC 



(4.9) 
(4.10) 



5 Gluon-parton collision 

We consider now gluon-parton scattering. In this case we have to distinguish two different 
final states for three-particle production. In addition to the parton b there may be: 

i) two gluons; 
a) a quark- antiquark pair. 



5.1 gb ggb 




Figure 7: Real gluon emission amplitude in gluon a fragmentation region. 

Assume that the additional gluon is emitted in the gluon a fragmentation region. The 
scattering amplitude in the Feynman gauge is given by 

Ay/ , _ ^/^i JM3 */^5 /taiagasfe J_ 9 1/ .b r 

1^2 
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where the amphtude A, corresponding to the diagrams of Fig. (H) was found in Ref. 
to be 



'P3 ^ P5 
+ I /is ^ /Us 

.03 ^ (35 , 



/il ^ /is 

ai ^ as. 



+ 



(5.2) 



where the current 



x-y 



P2-y\ ^ , P2-y^ \ u , t \ u 

P2-Z J X-Z 



(5.3) 



is a function of the momenta x", z", with x-y = x'^ya- 



Figure 8: Gluon to gluon fragmentation tensor A'Jj^^'^^^^^^. 




The squared hehcity amphtudes corresponding to Eq. (|57^ ) were given in Ref. 0. 
Since we work in D = A + 2e dimensions, we perform exphcitly the polarization sum 



on Eq. (5J-) in App. The averaged squared matrix element that we obtain is actually 
independent of the spacetime dimensionality (i.e. £- independent) so that we finally obtain 
the same formula as in Ref. ||TT| : 



dZid[ki]d[k2] ^ ^ ^ ^ ^ '^r(l-£)/i2- 



X 



CAas zlkj + (1 - zifq^ + 2:1(1 - zi)ki-q 
TT q^{z,k, + {l-z,)q? ■ ^'^ 

This expression explicitly exhibits the symmetry in the exchange — fci q ; Zi 1 — 

due to the identity of the gluons emitted in the fragmentation region of gluon a, but will 

be used for the softer gluon only {zi < |), in order to avoid double counting. 
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The function 

is related to the gluon to gluon sphtting function, and is factored out in Eq. ( |5.4| ) in the 
whole fragmentation region, as for quark scattering. In this case also the cross section 
matches the leading expression (|2.9| ) in the central region zi <^ 1, and thus is taken to be 
valid in the half phase space (|4.4|) in which two of the three gluons have positive rapidity. 



5.2 gb qqb 

The main contribution to the cross section from this kind of final states is reached when 
the fermion pair belongs to the same fragmentation vertex, as in Fig. (|5.2| .a). Other graphs 
like Fig. (|5.2| .b) in the case b = g are suppressed by a factor \t\/s because of fermionic 
exchange. 




Figure 9: Quark antiquark emission amplitude for b = g. a) The pair is emitted in quark a 
fragmentation region; h) This kind of diagram is suppressed by a factor \t\/s. 

Assuming qq being emitted in the fragmentation region of gluon a, and labelling q 
with "3" and q with "5" (ps = q), the corresponding amplitude is 

A^qqb = U3B;% 1 9 2pf^ tl 6x,x, (5.6) 
K2 

where is the sum of the diagrams depicted in Fig. (|IOD . 
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Figure 10: Gluon to quark fragmentation tensor B'^'l,. 



This calculation is very similar to the qb qgb cross section, and yields (App. 



do", 



qqb 



dzi d[ki\ d[k2 



in which the function 



a 



X 



klfi^e 2r(l - e)n^' 
CAas -2:1(1 - zi)ki-q 
n (ziki + (1 - Zi)q) 



TT 



(5.7) 



2Zi{l-Zi) 

l + e 



(5. 



related to the gluon to quark splitting function, is regular in zi G [0, 1]. Hence these final 
states do not produce any log s term, and need no subtraction. 

The virtual correction to the amplitude, including NLx effects, is given by an ex- 
pression analogue to Eq. (|4.5| ) where rit^ is replaced by the gluon-gluon-reggeon vertex 
correction in the helicity conserving channel 



r« = -«(fc?) 



1 11 

i ~ 12 



36 ~ T J iV^ U ~ 18^ 



(5.9) 



and provide the virtual contribution to the cross section 
-hf{k,)5{l-z,)5[q]2T^^h^\k,) 



(5.10) 



d^ld[fci]d[fe2] 

The total fragmentation vertex for the gluon is then obtained by summing Eqs. ( p. 41) , 
(|5.71 ) and (|5.10|) , dropping the common factor h^^\k2). 



In order to extract the gluon impact factor, we have to eliminate the logs term by 
subtracting the 1/zi singularity coming from the gluon splitting function ( ^.5] ) by means 
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of Eq. ( |3.4|) . Having identified in Sec. zi as the momentum fraction of the softer 
gluon in gb ggb scattering, to avoid double counting due to final gluons identity, we 
restrict the naive half phase space zi G [^, 1 — -^] to zi G [-^, |]. In so doing, the 1-loop 
correction to the gluon impact factor, explicitly computed in App. ^, is given by 

A 6 3ivJ + U m ) m V3 ^ 3V iVc J ^ ^ ^ 

where K, has been defined in Eq. ( |4.1CI| ). 



6 Collinear factorization and finite parts 



We want now to analyze the structure of the quark and gluon impact factors derived in 
Sees. and @. The explicit expressions reported in Eqs. ( ^^ ) and ( p. Ill) show the 
presence of 

• a /3-function coefficient —-h = 

• the finite part PbaiuJ = 0) in the u expansion of the Mellin transform of the a ^ b = 
q,g splitting functions for the incoming parton a under consideration, as depicted 
in Tab. (g); 



a common and e-finite factor 



This suggests writing the complete one-loop impact factors /la(fc) = hf\k) + h^p{k) for 
both a = q, g in the following manner: 



1-6 



a., / k 



e \/i 



,2\ £ 



2nbh^^\k) -J2h^c\k)P,, - h^^\k)lCe 

c 

hf{k)-{Y,hf\k)P^, + hf{k)lCe 



(6.1) 



+ 0{al) . 
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The factor in front of the last expression provides the renormahzation of the couphng 
constant, and by introducing the one- loop running coupling 



(0) 



2\ £■ 



(6.2) 



it can be incorporated in ha itself. 

Therefore, the impact factors in Eq. ( |6.1D assume the form 



as ^\ , as 



/,(o):^ + /,{o):^ 

^ e ^ e 



2\ e 



(6.3) 



and thus satisfy the DGLAP equations with splitting functions Pca{oJ = 0). 

After factorization of the collinear singularities contained in the latter, the only finite 
renormahzation is the one implied by the factor (^-'- + ^ ^) ' which is universal, i.e. 
independent of the parton type. 

Several comments are in order. First, the collinear singularities in do not contain 
the ~ terms of the splitting functions, which have been subtracted out in the leading 
term. In this context, we remark the difference between the double-fc cross section defined 
here through fc-factorization, and the double- minijet inclusive cross section defined by 



Mueller and Navelet [p!q] . 

In the latter case the cross section is inclusive over all fragmentation products of the 
incoming partons which are not identified, and measures the gluonic fc's only because 



a 



/o'z-Pca(^,£) dz 



q q 



Cf 



+ O{u) 



g q 



Cf 



1 + 



- ez 



g g 



2N, 



1- z 



z(l-z)--<5(l-z) 



2Nr 



nNc + {2 + e)nf 



q g 



2z(l - z) 
l + e 



2 1 



Oico) 



Table 1: Partonic splitting functions and Mellin transforms in uj space. 
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of collinear strong ordering. Therefore, it factorizes the gluon structure function in the 
parton, rather than the impact factor, with all its collinear singularities included and, 
furthermore, can be applied only at leading log s level because of strong ordering. 

On the other hand, in our case the definition of the fc-dependent cross section is more 
precise, but can be done theoretically rather than experimentally, because one needs not 
only to identify all fragmentation products, but also to subtract out the central region 
tail - which appears to be hard to do experimentally. Because of this difference, the 
collinear singularities to be factored out are different in the two cases. 

As a second point we remark that a universal renormalization, with the same /C coef- 
ficient, holds also for the soft part of the one-loop timelike splitting functions |2^. This 



part has a next-to-leading 1/(1 — zi) singularity, the leading one having a logarithmic 
factor, and is therefore analogous to the impact factor, which corresponds to the NLx 
constant piece in the high-energy limit. 



The above analogy is perhaps a hint [|15| to explain the universality found here. How- 



ever our result does not seem to follow in a clearcut way, because of the difference between 
collinear and high-energy factorization pointed out before. 



7 Discussion 

To sum up, we have defined here in a more explicit way the fe-factorization scheme intro- 



duced by one of us |[TT[], and we have applied it to the calculation of the finite one-loop 
corrections (Eq. ( |6.3| )) to the partonic impact factors, which turn out to be universal. 
Furthermore, we have generalized the definition of impact factors and of fc-factorization 



scheme to hard processes with colourless probes (Eqs. ( p.9|) and ( |3.1CI| )) in such a way as 



to be sure to avoid spurious infrared singularities in the subtraction of the leading terms. 

The above one-loop results, plus further use of fc-factorization and of the form ( |3.1D 
of the gluon Green's function, allows to extend the analysis of the NLx fc-factorization to 
two-loop level, and to extract the NLx BFKL kernel |1^, |13[, for a factorized scale of the 
type So = kik2. 
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In order to close the NLx calculations, the remaining task is to justify Eq. ( p. I]) to all 
orders. In fact, checking Eq. (|2.1|) at two-loops is not a real achievement, because it only 
determines the NLx kernel, depending on the scale choice. 

The question is then whether the kernel iterates correctly at higher orders for the class 
of factorized scales which is used. Investigating this problem probably requires the use of 
group-theoretical variables pO[ in order to incorporate exact s-channel phase space in the 
Mellin-transforms. It is hoped that this investigation will lead to a better understanding 
of the role of fe-dependent scales in the NLx kernel at higher orders as well. 



A Squared matrix elements in the fragmentation re- 
gion 

In this Appendix we report the explicit calculations of the gg* fragmentation vertices, 
both for gg and for qq production. Let's remind that pj, Xj and aj label respectively the 
momentum, the polarization and the colour quantum number of the j-th external particle. 

A.l Polarization vectors 

Here we present the conventions for the polarization vectors of the gluons that we have 
adopted in this paper. Even if it would be simpler to use real polarization vectors, we 
consider complex ones that in the D —>■ 4 limit reduce to the usual helicity basis. 

The light-cone coordinates are defined to be 

(i) - ( 'T^ -T*^ - -7*^^ I rf~- /v.^ /v*^ . \ 



' = 7= — ; X 



V2 

where x^ represents the extra coordinates in a 2e dimensional space. We introduce also 
the complex transverse coordinate x = x^ + ix"^. 

There are 2 + 2e physical polarization vectors and we take, in the p2 gauge, 

e(+)(g) =e(^)(g) = i=(0, 4; 1, 0) ; e(,)(g) = (0,0;0,0;r)) 
V 2 T 
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where 17 is the rj-th unit-vector in the 2£:-dimensional extra space: r/ G {1, ■ ■ ■ , 26}. Indi- 
cating with A any of the physical polarization, A G {+, — , 1, ■ ■ ■ , 2e} we have 

e(A)(g)-g = ; e(A)(g)-p2 = VA = +, -,77 . 

Defining the involution A 1— > A* acting on the polarization indices in the following way: 

(+r = (-) ; (-r = (+) ; {vr = {v) k^) = H^') > 

the orthonormality relations assume the form: 

e(A,)('?a)-e(Af,)(%) = S\aXb 

e(A.)(ga)-e(A,)(gb) = e(A,)(ga)-e(A,)(gfe) = -5a,a- = -5asa, yqa,qb ■ 

The polarization vectors {^(x){q) : A = +, — , 77} are a basis for the subspace orthogonal to 
the plane {p2, q) and 



A 

where n is the gauge vector {n = p2 for us). Contracting the Lorentz indices we obtain 



A. 2 gb ^ ggb squared matrix element 

The transition amplitude for gluon-parton scattering with production of an additional 
gluon in the fragmentation region of the incoming gluon a (Fig. (|^)) in the high energy 
limit s — > 00 has been given in Eq. in terms of the fragmentation tensor ( |5.2| ) 

'^^'fiTi^B = ^)ii%%i2p2 and the current (p^. By using physical polarization e-p2 = we 
can drop the p2"Component of D^^ and retain 



1 



x-y 



P2-y 



D^'{x,y,z) = (y-z - P2-P4)x'' + {x-p2)y 



x-z 



Applying in Eq. ( |5l^ ) the Jacobi identity 
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the fragmentation tensor reads 

- ^A.iM3^M5(P3, -PuPb) - ^M3M5^Mi(P3,P5, -Pl)] + (3 ^ 5)| 

where, for brevity, we have defined 

D^ = D^(-P1,P5,P3)+^'^(P5,-P1,P3) , 
D^=D^(P3,-P1,P5) , 

D^c = D''{P3,P5,-Pi) . 

By averaging over all colours and polarizations, we obtain the unpolarized emission prob- 
ability 



^ggb - 



1 1 



^3 ^5 ^, 



aiaaas b 



(iv2 - 1)2 2 + 2£ my 

Now we isolate colour and polarization structures by defining 

^135 = Zx,X;X,{Pi,P3,P5) = e^*^ ef {gf,,^,,D^, + ^mims^Js + ^msms^S) 

in such a way that 



and 



(2 + 2e) (iV,^-l)^ 

16^6 CbiV,' 1 
2 + 2e - 1 {kjy 



}^ \ ^ l^aiase^easb ^ _|_ ^aiase^easb ^ I 

.2\2 / ^ I I 



|^135p + |-Z'l53p + o(-^135-^153 + -^153-^: 



135 j 



(A.2) 



We evaluate the terms in square brackets by using Eq. ( |A.1D . In this way we get 



J2 1^1351' = - (2 + 2e) \Ao,p{ps,P3)D'XD^^ + A^p{p,,p5)D%D^^ + A,^(pi,pi)DgDg 



+ 2 



A^p{p,,p,)D'XD^^ + A^p{p,,p^)D'XD'^ - A^p{p,,p,)D^,D 



a r,/3 



B 



(A.3) 
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By expanding the Lorentz contractions substituting the appropriate kinematical variables, 
the interference terms (i.e. the second hne of Eq. ( |A.3| )) cancel each other, and the first 
line yields 

The second term in ( |A.2|) has the same structure as the first one, and gives 

21 

'Pi^Ps I fci 

Aj Aj \j 

For the third term in (|A.2| ) we obtain 



(A.5) 



X ^1*35^153 =(2 + 25) [a„^(p3,P3)I^A^B + A,^(p5, P5)I^S/^1 " A„^(pi, 



Also in this case the interference terms cancel each other, and we obtain 
2^^135^53- Z:2-J7-7-— TT^-T-Ti • 



The fourth and last term in (|A.2|) is the complex conjugate of the third one and hence, 
since ([A.6| ) is real, they are equal. 

Multiplying Eq. ( |A.2| ) by the three-body phase space ( p.8|) and substituting the ex- 
pressions ( |A.4| , [A.5| , |A.6|) , we finally obtain the differential cross section ( |5.4| ) for gb — > ggb 



scattering. 



A. 3 gb qqb squared matrix element 

The calculation of the gb qqb cross section is analogous to that of qb ^ qgb and 
requires the evaluation of the gg* qq vertex, which can be obtained by means of 
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fc-factorization 



(2 + 2^)(^e^ 
Adopting the Feynman gauge, 

and by using physical polarization with e-p2 = and the high energy kinematics, we can 
write the colour decomposition 



where 



Pig 

,(-) 12 1 2 
fel ^ = + 62 = + • 

PiQ Psq P1P3 Psq 

By performing the gamma matrices algebra and the polarization sum in Z) = 4 + 2^ 
dimensions, we obtain, for each colour structure, 

J2 lusBv^l^ =(2 + 2e)s^y^g{zi,e) x 
A.- 

X {bj{l - zi)piq + bjzipips + hih2[pzq - ziPiPa - (1 - 

where CPqg is given in Eq. (^.8|) and is related to the gluon to quark splitting function. 
By averaging the square of the amplitude ( ^.61) over all colours and polarizations, the 
emission probability reads 

TT 1 Cb 



{klY - 1 

which, multiplied by the three-body phase space ( |2.8| ), gives the differential cross section 
( pT7| ) for gb qqb scattering. 
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B Calculation of ki integrals 



In this Appendix we show the relevant integrations (over zi and in ki transverse space) 
of the fragmentation vertex in order to get the 1-loop impact factors, both for the quark 
and for the gluon case. 



B.l Quark impact factor 

In first place, we show that the Sudakov term of qgb production, upon ki and zi integra- 
tion, cancels out with the Cp term of the virtual correction. In fact, the Cp contribution 
to the fragmentation vertex Fq{zi, ki, ^2) is 
?gq(^i,e) Cpas zfkl 



r(l-£)/x2e ^ q^{q-zik2 



2 

TT 2Cf01sN^ CpOis 



iV2-l //25 27rr(l - - L 'v ' '^q^{q-Zik2Y 

In order to perform the fci-integration, we use the formula 

y r(a + /? + 2£)r(i-a)r(i-/?) ' 



to obtain 



/ d[fci] f d[g] _ V{l-e)T\e) 

J q'iq-z.k^y J q'iq-z.k^y r{2e) ^ 



2\e-l 2e-2 
2) ^1 



The integral over zi, due to the vanishing of the integrand when zi goes to 0, can be 
extended in its lower limit from q/ a/s - the half phase space - down to 0, introducing in 
this way a negligible error in the limit s — > 00. This yields 

7V2-1 27r//2^^''2j p(2£) io ' ^ ^ ^1) 



Cpa, T\l + e) fklVfl 2 1 



^ TT eT{l + 2e) J \e l + 2e ' 2 
= /.('^)(fe2)2a;W(fe^)^(-i + ^-4e)+0(£) 

which exactly cancels the part of the virtual term proportional to C^r 

d[k,] fdz, h^^\k,)2T^^\ 5{l-zMQ] = hf{k2)2uj^\kl) rW 
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as can be seen with a glance to Eq. (|4 



The effective quark fragmentation vertex is then as given in Eq. ( [4.7[ ), i.e. 



F^izu k^, k2) = h^^\k,) 5[q]5{l - z,)2cu^'\kl 



11 



12 



85 TT^ 



36 



iVe V 6 1^ 



Tgq(2i,e) as (l - zi)q-{q - Zik2) 



[q - Zik2Y 



r(i - £)/i2^ g2 

In order to obtain the 1-loop correction to the quark impact factor, we have to use Eq. 
to subtract the 1/ Zi singular part, which gives the logs behaviour of the cross section. In 
this way we get 



a 



,F^'\z,,k,,k2 



2CFasNe X 



(B.2) 



%,(1-.)^ 



11 

12 



85 vr^ 
36 4 



e — 



N,\Q li 



a. 



2z^'^^' ^'''^ klqHq-z,kl) z.k^ 



r(i-£)/i2^ 

and then we have to integrate over zi G [0, 1] and ki at fixed ^2. 

The integration of the first term (i.e. the virtual contribution) is straightforward and 
yields 

'85 7r2 
18 ^ 2 



hf{k2)uj^^\k 



11 



Uf ( \ 5 

' e 

3 9 



(B.3) 



Using the identity Q{zi) = 1 — 0(— 2:1) the last line of Eq. ( p.2|) yields the multiple integral 



a. 



ni-e) 



dzi 



?gq(Zl,£)(l-Zl) / d[k 

/ki 1 
d[ki] log— Qk^q -72- 
q k^q 



q {q - Zik2) 



1 



d[fci] 



klq'{q-zik2f zj klq\ 

(B.4) 



It is convenient to perform in first place the ki integration, since dimensional regulariza- 
tion provides simple integrals in the transverse momentum space 



The integrand of the first fci-integral in ( |B.4|) can be written as a sum of terms with 
two factors in the denominator: 



q {q - Zik2) _ 1 

klq^ {q- Zik2y klq^{q-zik2y 



^q' 



kl + 



1 



2(1 -zi) 

-v2 



2(1 -^i) 
1 



{q - zik2f 



2(1 - z^) [kl {q - z^k2Y q" {q - ^1^2)' kl q\ 
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thus allowing, with an appropriate shift in the integration variable, to apply (|B.1|) with 
a = P = and to obtain 

J ^ ^^Iq^q-z^k^Y eV{l + 2e) ^ ^' I - 

The second fci-integral is immediately solved: 

f d[fci] _ T{l-e)V\l + e) , , 
J k\q^- eT{l + 2e) '^"'^ 

For the third integral we use the representation 

, , /•+^°° dA 1 /a\A 

log - Bab = hm 



b a^o+ J_ioo 27rz {X + ay V6. 

valid for a,b > 0, which allows us to write 

/ d[k,] log ^ e.,, -r^ = ^ lim ^ / d[k,] (^) ' -1^ 



2a^o+y_,^ 27rz(A + a)2 r(2£)r(l - A)r(l + A) 



The function 



A ^ n>^ + e)n-X^e) 

r(2£)r(l-A)r(l + A)(A + a)2 ^ ^ 

rapidly vanishes for |A| — > oo in all directions, apart for the two real semi-axes. Therefore, 
we can displace the integration contour aroun the positive real semi-axis, enclosing the 
poles placed inA = e + ?T, inGN. The residue of ( [B.5D in A = e + n is 

(-1)"+^ T{n + 2e) 

n\ T{l + n + e)V{l-n-e){n + a + ey 

In this way we can set a = and replace the integral with the sum 
ir(l-5)„2,e-i^(-ir T{n + 2e) 



2 V{2e) ' ^0 ^- r(l + n + e)r(l-n-£)(n + e)2 

The first term in the sum (n = 0) is 

r(25) 
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which is 0{6 ^) and contribute to the integral to 0{6 ^). The other terms (n ^ 1) are 



(^(e^) and contribute to the integral to 0{e). In conclusion 



T{l-e)T\l + e) 1 



_ (k'^y-^ 

q -■ kiq^ eT{l + 2e) 2^ ^ 



0{e) 



(B.6) 



We are left with the zi integration: 



eT{l + 2e) ^"^^ 



a 



X 



1 

2x1 



l + {l-z,f + ezi][l-zl' + {l-z^) 



\2e 



2 

Zl 



1 

2l 



We observe that the 1/2:1 singularity coming from ^Pgg is correctly subtacted out, and the 
l/2e term produced by the log(A;i/g) integration cancels out with 

2 



ziY + ezi 



_ + ^(l)_^(l + 2e)-l+ ^ 



1 - zf + (1 - z^f' 
l + e fl 1 



1 



Summing the virtual contribution 



2 l + 2e l + e^ 
and 2CFasN^fi~'^'^Iq yields finally 



a 



h^\k,)uj^^\kl 



"q K'-s)- \-''2) 

as given in Eq. (|4.9|) of the text. 



11 



3 1 
2-2^ 



67 

18 ~ ¥ 



5^/ 



B.2 Gluon impact factor 



The subtracted fragmentation vertex F^'^ for the gluon, calculated in Sec. (H), is 

~ 2CAasN, 



asF^'^\zi,ki,k2 



X {6[q]6{l- Zl) 



+ nf 
+ 



1 11 /67 7r2\ Uf /I 5 
e ~ 12 ^ V36 ~ Ty ^ ^ iVc V6 ~ 18^ 
-2:1(1 - zi)fcrqr _^ Cira^ 



yqg(^i,g) [Os 
4r(l-e)/x2s [tt klq^Zlkl + {l- Zi)qy ' CaTT fc^ g2 

a.. 



2:2fc^ + (1 - 2:i)2q2 + ^^(1 _ 



r(l-e)/i2^ 2ifc2g2 
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Let's now integrate over zi and ki in order to obtain the one-loop correction to the gluon 
impact factor. 



The first hne in curly brackets, corresponding to the virtual contribution, yields im- 
mediately 



'2 11 /67 TT^\ Uf /I 5 



(B.7) 



The second line corresponds to the qqb final state. Since it is finite in the limit 
Zl 0, we may extend the Zi integration from the half phase space Zi G 1] down to 
Zl = introducing in this way a negligible error in the s ^ oo limit. To perform the ki 
integration, we decompose the bigger fraction in a sum containing terms which have only 
two factors in the denominator, allowing us to apply Eq. ([B.l|) and to obtain 



h'f\k2)uj^'\k 



n,/ 1^23^ 



Nr 



nfCp (2 1 

m V3 3 



(B.8) 



The third and the fourth line represent the ggb contribution and the central region 
subtraction respectively. Let's define, for brevity, the function 
zfkl + (1 - Zi)2qr2 + ziil - Zi)ki-q 



R{z.ki,q)^ 
Using the symmetry properties 



;i - zi)q) 



R{1 - zi,ki,q) = R{zi, q, ki) = R{zi, -q, -ki) 
and Eq. i f>.5\) which allows us to write 

'?,,izi) = Pizi) + F{l-Zi) ■ P(^0=(^ + |)(l-^i) , 



we have: 



dzi d[ki]'>e{- ~ ZiP,,{zi)R{zi,ki,q) 



9(g - ziki 
k\q' 



Zl 



f dzijd[ki]!^P{zi)R{zi,ki,q) 



+ 72-^ log — ©fci^ 



klq^ Zl k\q^ " q 

where we have changed variables Zi ^ 1 — Zi and —ki fci + ^2 in the intermediate 
steps. The function P(-2i) has only one pole at = which is correctly subtracted by 
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the term. Decomposing R{zi, ki, q) in a sum of terms with only two factors in the 
denominator, and performing the ki integrations as explained for the quark impact factor, 
we get, at 0{e), 







dzi /d[fci] e(- - z^)?^,{z^)R{z^,k^,q) 



6(g - ziki] 



Zl 



2\ 



TT 2CAa,N, as r\l + e) ^ k^ ^ i 



2 11 /67 2n^ 



h^'\k,)co^'\kl) 



p{z,){i+zr+{i-z,r)--+^ 

Zl 

(B.9) 



Summing Eqs. ( |B.7| ), (|B.8|) and (|B.9|) we finally obtain the one- loop correction to the 
gluon impact factor: 



h<i\k2)= j A[ki] j\ziF^^\ziMM 



11 _ J!^ V ^ii + + 
6 ?,NJ V 6 QN^ 



3 + 



as given in Eq. (|5.11|) of the text 
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